Controlling diffusive transport in confined geometries 
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We analyze the diffusive transport of Brownian particles in narrow 
channels with periodically varying cross-section. The geometrical confine- 
ments lead to entropic barriers, the particle has to overcome in order to 
proceed in transport direction. The transport characteristics exhibit pecu- 
liar behaviors which are in contrast to what is observed for the transport 
in potentials with purely energetic barriers. By adjusting the geometric 
parameters of the channel one can effectively tune the transport and dif- 
fusion properties. A prominent example is the maximized enhancement of 
diffusion for particular channel parameters. The understanding of the role 
of channel-shape provides the possibility for a design of stylized channels 
wherein the quality of the transport can be efficiently optimized. 

PACS numbers: 05.60.Cd, 05.40.Jc, 02.50.Ey 

1. Introduction 

Diffusion process in narrow confined systems exhibits peculiar proper- 
ties which are radically different from what we generally observe in free 
systems, i.e., in the absence of any geometrical restrictions [I]. If we con- 
sider a Brownian particle which is moving in a potential energy landscape 
or in a confined geometry where the geometrical confinement may regu- 
late or control the diffusion process, the diffusion varies significantly from 
the free case [21 OH]. Depending on the conditions imposed, the diffusion 
coefficient can be larger or smaller than the bulk diffusion coefficient. Brow- 
nian particles when moving in a confined geometry undergo a constrained 
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diffusive motion. This feature of constrained motion is ubiquitous in ion 
channels, nanopores, zeolites, and for processes occurring at sub-cellular 
level [H El El [3 El El [TOl HH [ISl H3l HH HSl US]. The uneven shape of these 
structures regulates the transport of particles yielding important effects ex- 
hibiting peculiar properties [41 113 1 1^4" ! I15j. The results have implications in 
processes such as catalysis, osmosis and particle separation [51 \6[ [71 [8| \TT\ [T2] 
and on the noise-induced transport in periodic potential landscapes that lack 
reflection symmetry, such as ratchet systems [TTJ [THJ Q32 [2D] . 

Nowadays, artificial, i.e., synthetic, nanopores are available [3 EH [22] to 
characterize the transport characteristics of ionic species like K + ,Na + , CP 
[71 HH] and of macro molecules, like for example, DNA or RNA [221 [231 IM [251 
[26] . With latest technology, these synthetic nanopores can be made upon 
choice, allowing for effective control of the diameter and the shape of the 
nanoporous structure [71 l21[ [22]. In recent years, it has been of great interest 
to reveal the sequence and structural analysis of DNA and RNA molecules 
by passing them through nanopores [23j [Ml US [26J [271 [28]. When a double 
stranded DNA molecule passes through the charged nanopore, each base 
pair exhibits its own distinct electronic signature because each base pair is 
structurally and chemically different. While present conventional methods 
to find out this sequence would take several months and are expensive, 
one could find the sequence of a human genome in a matter of hours at a 
potentially low cost using this new technology. In addition, charging of the 
inner tube walls could lead to rectification of the motion of ions [TJ [161 EH] • 

It has been found that the separation of macromolecule fragments and 
ions in narrow channels [1311141 IT5], I22 |, I24j is largely influenced by their shape 
alone. The characteristic behavior of ionic currents of these molecules or 
individual species can effectively change if the pore diameter changes [Ml 
|2"B"1 [3U1 EH [32]. The subject of this paper is to systematically illustrate the 
role of the channel properties on the transport characteristics for Brownian 
particles passing through the channels. 

The paper is organized in the following way. In section II we introduce 
the model system, the dynamics and discuss the complications involved in 
solving the full problem. The reduction of dimensionality, or in other words 
the simplification, and the analytical treatment of the problem is discussed 
in Section III. In Section IV we introduce the simulation techniques to 
solve the full problem numerically. The main findings on the transport 
characteristics are discussed in Section V. In Section VI we present the 
main conclusions. 
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2. Modeling 

We investigate the transport of Brownian particles through straight 
channels with periodically varying cross-section. In the following, we mainly 
focus on 2D - channels (see Fig. [1]) , although a very similar line of reasoning 
could be applied for pore structures in 3D. For a model system we consider 
a channel structure whose shape is defined by the periodic boundary func- 
tion uj(x) with the periodicity L, i.e. oj{x + L) = uj{x). The diffusive motion 
of the Brownian particle is then confined by the upper boundary u(x) and 
the lower boundary — lo(x). In particular, we choose: 

u(x) =^sin(2vry) + J B, (1) 

Li 

with the two channel parameters A and B and where the function u(x) 
is thought to be the first terms of the Fourier series of a more complex 
boundary function. Due to the symmetry with respect to the x-axis the 
boundary function could be given in terms of the maximum half-width of 
the channel w majX = (A+B) and the aspect ratio of maximum and minimum 
width e = (B — A) j{B + A), i.e. 

W (x) = ^i(l- e ){sin(2^) + i±l}. (2) 

Transport through the considered confined geometry may be caused by 
different particle concentrations maintained at the ends of the channel, or 
by the application of external forces acting on the particles. Here we will 
exclusively consider the case of force driven transport. The external driving 
force is denoted by F = Fe x which points into the direction of the channel 
axis. Then the dynamics of a suspended Brownian particle is described by 
means of Langevin equation, in the over damped limit [33] > which reads 

H T -» , — * 

7 _ = F + v /^r£(t), (3) 

where t is time, x the position vector of the particle, 7 its friction coefficient, 
&b the Boltzmann constant and T the temperature. The thermal fluctuat- 
ing forces which the surrounding fluid exerts on the particle are modeled 
by zero-mean Gaussian white noise obeying the fluctuation-dissipation 
relation {£,i(t) £,j(t')) = 25ijS(t — t') for i,j = x,y. In addition to Eq. ©, 
the full problem is set up by imposing reflecting boundary conditions at the 
channel walls. For the sake of simplicity, we measure lengths by the period- 
icity of the channel L, time by r = 7L 2 / (/cbTrJ which is the corresponding 
characteristic diffusion time at an arbitrary but irrelevant reference temper- 
ature Tr, and force by Fr = /cbTr/L. Consequently, the Langevin equation 



Fig. 1. Schematic representation of a single channel with the periodicity L, bottle- 
neck half- width w m i n and maximal half-width ui max . The shape of the structure is 
defined by Eq. @. The constant applied bias F acting only along the length of 
the channel. 

reads in dimensionless form: 

^ = Fe x + y/Dg(t), (4) 

where D = T/Tr is the dimensionless temperature, and the boundary func- 
tion is given by: 

W (x) = (1 - e) {sin(2™) + I±f } . (5) 

The corresponding Fokker-Planck equation for the time evolution of the 
probability distribution P(f, t) takes in dimensionless units the form |34}I35|: 

= -V ■ (F - DV) P(r,t) . (6) 

Since, in the present situation, we deal with irregular and impenetratable 
channel walls, the reflection of particles at the boundaries leads to a vanish- 
ing probability current [1J . Solving the full problem in a higher dimensional 
space with this irregular boundaries is a difficult task and so far there is no 
analytical method. Despite the inherent complexity of this problem an ap- 
proximate solution can be found by introducing an effective one-dimensional 
kinetic description where geometric constraints and bottlenecks are consid- 
ered as entropic barriers [H [361 E2 EH1 EH] . 

3. The origin of entropic barriers 

The dynamics of the system can be approximatively described by means 
of a ID kinetic equation, obtained originally from 3D (or 2D) Smoluchowski 
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equation by assuming fast equilibration in orthogonal channel direction. 
This so-called Fick-Jacobs equation reads in presence of an applied bias [1] 



dt dx \ dx D dx 

and approximate the full dynamics for u(x) <C 1. In Eq. fl7|), P(x,t) is 
the probability distribution function along the length of the 3D tube or 2D 
channel, A(x) defines the free energy: 

A{x) -=E-TS = -Fx -Din s(x) , (8) 

where s(x) is the dimensionless transverse cross section s(x) := t:uj 2 {x) in 
three-dimensions, and dimension less width s(x) := 2ui{x) in two-dimensions, 
where oj(x) is the radius of the pore in 3D (or the half- width of the chan- 
nel in 2D). Interestingly, bottlenecks in the channel structure accounts for 
barriers in the potential function A{x) whose height scales with the tem- 
perature indicating a clear entropic contribution. In this terms our model 
system allows for investigation of transport in presence of entropic barriers. 
The spatially diffusion coefficient appearing in Eq. (J7|) reads 

where a = 1/3,1/2 for two and three dimensions, respectively. Here, in- 
troducing the spatially diffusion coefficient, into the kinetic equation, im- 
proves the accuracy, and extends it's validity to more winding structures 

gl E2 Eg E21 EQ1I1I1 021 ■ 

4. Transport characteristics 

For the transport of Brownian particles the quantities of interest are: the 
mean particle velocity (x), respectively the nonlinear mobility \i = {x)/F, 
the effective diffusion coefficient D e g, and a measure for the quality of the 
transport. For the latter we take the so-called Q-factor, which is the ratio 
of the effective diffusion coefficient and the particle current, reading 

«=W" (10) 

4-1. Analytics 

After applying the equilibration assumption pointed out in Sec. El Eq. [7] 
could be derived. For periodic structures analytical expressions for the 
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above mentioned transport characteristics can be determined by using the 
mean first passage time approach [H \39\ HO]. Accordingly, the average par- 
ticle current is given by 

1 - e~ F ' D 

(i) = — , (11) 

I(z) dz 



[I{z)\ 2 dxdz 

, (12) 





and the effective diffusion coefficient by 

r 1 r x D(z)e A ^/ D 

D eS = D 



o J x -i D(x)e A (*)/D 



i i 3 
I(z) dz 



o 

where the integral function reads 

,A(x)/D r x 



D(x) 



jr W = £ r5V /" dye- A{y)/D - (is) 
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In the case of an energy barrier, the driving force F and the temperature D 
are two independent variables, whereas for entropic transport, both current 
and effective diffusion are controlled by a universal scaling parameter F/D 
[B IH [39j SO] . As, for a given structure, the nonlinear mobility fi = (x)/F, 
the ratio of effective diffusion D e g and bulk diffusion D solely depend on 
the scaling parameter F/D, the same applies for the Q-factor and Eq. (|10p 
reduces to: 

Q = £^lR Qheei (14) 

where Qf ree is the Q-factor for the biased Brownian motion in absence geo- 
metrical constraints (free case) which is the inverse of the scaling parameter, 
i.e., Q free = 1/(F/D). 



4-2. Brownian dynamics simulations 

The approximative transport characteristics described above can also be 
compared with those obtained from precise numerical simulations consid- 
ering the full 2D dynamics. The mean particle current and the effective 
diffusion coefficient, have been corroborated by performing Brownian dy- 
namic simulations by integrating the Langevin equation Eq. ([3]), within the 
stochastic Euler-algorithm. Then, the mean velocity in x-direction is given 

by 

(x) = Urn ^ , (15) 



7 



and the corresponding effective diffusion coefficient reads 

i^-lh. <*■<'» -<*»' ■ (16) 

t— >oo Zt 

By comparison of the precise numerical and the analytical results for 
different channel geometries, the validity of the equilibration assumption 
could be analyzed in detail [39\ BP], 



5. Geometry controlled transport 

The confinement by the considered channel geometry can be altered 
by systematically changing the parameters w max and e in the geometric 
function, Eq. [5j Changing these parameters we can consider two cases: 
constant-ratio- scaling and constant-width- scaling. In the former case, we 
modify the maximum width w max of the geometry and keep the aspect 
ratio e constant. Consequently, by applying the constant-ratio-scaling the 
minimum width, i.e. the width at the bottlenecks, is changing by modifying 
the maximum width. In contrast, in the constant-width-scaling, we keep 
the maximum width w mSiX constant and change the aspect ratio e. In the 
following we analyze the transport characteristics within the two different 
scaling regimes. 



5.1. Constant-ratio-scaling 

In constant-ratio-scaling, we fix the ratio of maximum and minimum 
width. Hence, by increasing the in sxiniiin i width if max the channel is diluted 
in orthogonal channel direction. The advantage of this scaling is related to 
the effect, that within the approximative reduced dynamics, the height of the 
entropic barriers is kept fixed while the curvature of the effective potential 
at the minima and maxima varies. In Fig. [2] the mean particle current 
and the effective diffusion coefficient are depicted for different maximum 
channel widths. Strikingly, the particle current (x) decreases with increasing 
effective temperature D [H 0]. This is in contrast to the transport over 
energetic potential barriers where thermal activation leads to an increasing 
particle current when the temperature is increased [34J. In our system, 
temperature dictates the height of the entropic barriers. With increasing 
temperature, the height of the barriers increases and the particle current 
decreases, cf. Fig. [2|a). In the deterministic limit, i.e., D — > the average 
particle current equals to the applied bias F irrespective of the maximum 
channel width. In this limit, the particles do not explore the side bags of 
the channel structure and move within a region defined by the width at the 
bottleneck. Consequently, the particles assume the velocity of a biased free 
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Fig. 2. (Color online) Noise (temperature) dependence of the average particle cur- 
rent (a) and the effective diffusion coefficient (b) for a symmetric two-dimensional 
channel with the shape defined by Eq. ([5]) at the force value F = 0.628 and for 
various values of the maximum half-width (w ma x) of the geometry and for a con- 
stant aspect ratio e = 0.01. The lines correspond to the analytic, approximative 
results given by Eqs. (fTTjl and (TT2"]) . The different symbols correspond to simulation 
results, cf. Eqs. (fT5"j) and ([TBI) . The arrow indicates the particle current for the 
deterministic limit for which the particle current equals F. 

particle. However, this observation is not captured by the analytics, as the 
influence of the winding of the structure is overestimated in the spatially 
dependent diffusion coefficient [3D] . 

As the maximum width increases, the area of the size bags increases. 
More available space in the orthogonal channel, however, leads to reduced 
particle current along the channel direction. Therefore, the particle current 
decreases with increasing maximum width, cf. Fig. [2(a). A similar reasoning 
applies for the increased enhancement of diffusion, which becomes visible in 
Fig. E(b). The effective diffusion coefficient does not depend monotonously 
on the temperature D, and exhibits the effect of enhancement of diffusion 

EQ- 

Overall, for small values of w max , i.e., smooth geometries the analytical 
description leads to better results whereas it fails for large u> max and small D 
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Fig. 3. (Color online) Noise (temperature) dependence of the average particle cur- 
rent (a) and the effective diffusion coefficient (b) for a symmetric two-dimensional 
channel with the shape defined by Eq. §5§ at the force value F = 0.628 and for 
various values of the aspect ratio (e) of the geometry and for a constant maximum 
half- width u> max — 1 . Like in Fig. [2] lines correspond to analytic results and sym- 
bols to numerical results. Moreover, the arrow shows the particle current for the 
deterministic limit. 

range [U HQl [39] . Also the applicability of the analytical description depends 
on the parameter we look at, i.e., whether it is the average particle current 
(1st order) or the effective diffusion coefficient (2nd order). 

5.2. Constant-width- scaling 

If we keep the maximum width constant and increase the e-value, the 
strength of the entropic nature of the system, i.e. the height of the entropic 
barrier in our system, is reduced. In Fig. [3] we depict the dependence of the 
average particle current and the effective diffusion coefficient on the noise 
strength. Once more, the particle current exhibits this peculiar dependence 
on the noise level which is characteristic for the transport over entropic 
barriers. Namely, the particle current decreases with increasing noise level. 
In the deterministic limit, the particle current tends to the applied bias F 
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for all values of e, cf. Fig. E^a). Upon decreasing the aspect ratio e, the 
width of the bottleneck decreases and the effective height of the barrier 
increases. Consequently, the particle current decreases with decreasing e. 

Moreover, by decreasing the barrier height, i.e. increasing the e-value, 
the effect of the enhancement of diffusion becomes less significant, and in 
the limit of a flat channel, i.e. for e — > 0, there is linear, monotonic behavior 
of the effective diffusion coefficient, cf. Fig. [31(b) . 

Overall, for the constant-width-scaling, the analytical description leads 
to better results for the geometries with very small aspect ratio. With 
aspect ratio e, the bending of the boundary function and, therefore, its 
derivative increases and the criteria for the applicability of the Fick-Jacobs 
approximation, i.e. co'(x) <C 1, increasingly fails, cf. Fig. [3j 



5.3. Comparison of the different scalings 

The transport process depends on many factors such as the slope of the 
structure, the width of the channel at the bottleneck, the strength of the 
applied bias, and the thermal noise present in the system. The quality of 
the transport can be measured by the Q-factor which is introduced in Sec. HI 
A large Q-factor means more randomness in the transport process. 

Fig. H] shows the dependence of Q for the two different scalings. In- 
terestingly, with increasing maximum width of the channel, while keeping 
the aspect ratio of channel constant, the Q-factor increases. Means, the 
transport becomes more noisy for larger channel widths, cf. the behavior of 
Q-factor in the inset of Fig. Eta) . The maximum of the Q-factor for wider 
channels is accompanied by a increased enhancement of diffusion on one 
side and on the other side the reduced particle current, cf. Sec. 15.11 and, in 
particular, Fig. O 

Within the constant-width-scaling, smaller e leads to a decrease of the 
height of the entropic barrier. Hence, the transport should be better in the 
sense of a smaller Q-factor. However, we found that for rather small and 
fixed values of u> max an increase of e could spoil the quality of transport by 
increasing Q, cf. the behavior of the Q-factor in the inset of Fig.^b). Thus, 
by varying the aspect ratio e of the geometry one can encounter different 
transport regimes. If the particular scaling parameter F/D is given, one 
can tune the quality of transport by choosing an optimal geometry, i.e., an 
optimal e-value. On the other hand, if the structure is given, the scaling 
parameter F/D may be used to control the transport. An optimal transport 
in the sense of a smallest Q-factor could always be achieved at either very 
small or very large scaling parameters F/D, cf. Fig. EJ 
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Fig. 4. (Color online) The behavior of Q factor as a function of the scaling pa- 
rameter F/D obtained from numerical simulations (a) for the case of constant- 
ratio-scaling, where the aspect ratio e = 0.01, and for various values of the max- 
imal half-width ui max . Circles: w ma x = 2.02; triangles: w max = 2.02/2; squares: 
«) nial = 2.02/27T. The inset in (a) depicts the behavior of Q as a function of w max 
at F/D = 10. (b) for the case of constant- width-scaling, where the maximal half- 
width w max = 1, and for various values of the aspect ratio e. Circles: e = 0.001; 
triangles: e = 0.075; squares: e = 0.03. The inset in (b) depicts the behavior of Q 
as a function of e at F/D = 10. The dashed line in both the figures correspond to 

Q free - 1/{F/D). 



6. Conclusions 

In this work we analyzed the transport of Brownian particles through 
channels with periodically varying width, which exhibit narrow bottlenecks. 
The effect of confinement can be recast in terms of an entropic potential, 
and the dynamics of the system can be effectively described by means of 
the Fick-Jacobs equation with a spatially dependent diffusion coefficient. 
We observe peculiar transport phenomena, which are characteristic for the 
transport over the entropic potential barriers. The transport characteristics 
like the average particle current and the effective diffusion coefficient sensi- 
tively depend on the geometry of the channel. Depending on the geometric 
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parameters we have considered two scaling regimes: the constant-width- 
scaling, for which we keep the maximum width constant, and the constant- 
ratio-scaling where the ratio of the maximum and the minimum width of the 
channel is kept fixed. For the former case, the diffusion process varies upon 
changing the aspect ratio of the geometry. There is a critical bottleneck 
opening related to a maximum effective diffusion coefficient. This feature 
is more visible in highly confined geometries, i.e., for smaller w m3jX values. 
In both geometric scalings, during the crossover region from a purely en- 
ergetic to entropy dominated regime, optimal transport could be observed, 
suggesting that by increasing the noise strength D the system may exhibit 
better transport, i.e., smaller Q— factor. Thus, tuning the parameters the 
aspect ratio e, the maximal half-width w mSLX , and D one can effectively reg- 
ulate the transport characteristics in these confined geometries. However, 
in general, controlling the width ratio 6 and the maximal half-width iy max of 
the geometry may not be feasible, but, tuning the scaling parameter F/D 
one can arrive at an optimal transport regime. 
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